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It has been shown that many of the phenomena related to the formation of "tails" in the low- 
concentration region of ion-implanted impurity distribution are due to the anomalous diffusion of 
nonequilibrium impurity interstitials. These phenomena include boron implantation in preamor- 
phized silicon, a "hot" implantation of indium ions, annealing of ion-implanted layers et cetera. 
In particular, to verify this microscopic mechanism, a simulation of boron redistribution during 
low-temperature annealing of ion-implanted layers has been carried out under different conditions 
of transient enhanced diffusion suppression. Due to the good agreement with the experimental 
data, the values of the average migration length of nonequilibrium impurity interstitials have been 
obtained. It has been shown that for boron implanted into a silicon layer preamorphized by germa- 
nium ions the average migration length of impurity interstitials at the annealing temperature of 800 
°C can be reduced from II nm to approximately 6 nm due to additional implantation of nitrogen. 
The further shortening of the average migration length is observed if the processing temperature is 
reduced to 750 °C. It is also found that for implantation of BF2 ions into silicon crystal, the value 
of the average migration length of boron interstitials is equal to 7.2 nm for thermal treatment at a 
temperature of 800 °C. 



PACS numbers: 61.72.Ji, 61.72.Tt, 61.72.Ss, 85.40.Ry 
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I. INTRODUCTION 



It is well known that boron is a basic impurity of p- 
conductivity used in the technology of the production of 
silicon integrated microcircuitsi. Unfortunately, boron 
atoms have a small mass and large mobility in silicon 
crystals. Due to the small masses of boron ions, the 
formation of an amorphous layer in the ion-implanted 
substrates is not observed even at large fluencies. Only 
a great number of radiation defects are created. The 
absence of an amorphous layer and the presence of radi- 
ation defects result in the significant transient enhanced 
diffusion (TED) of ion-implanted boron during the sub- 
sequent annealing (see, for example^'"— ). All these 
phenomena substantially complicate the problem of the 
formation of very shallow junctions with high electro- 
physical parameters. For suppressing the TED of ion- 
implanted boron, a method of boron implantation in a 
silicon layer preamorphized by heavier germanium ions 
is widely used^ii^"— !^"— . Due to the solid phase epitax- 
ial regrowth (SPER) of the amorphous layer, the region 
doped with boron is characterized by a perfect crystal 
structure, containing defects that are invisible by electron 
microscopy. However, the transient enhanced diffusion is 
observed as before, although it has another character and 
a smaller intensity. In the papers'^^''^'^ a qualitative differ- 
ence of the form of boron profiles produced by annealing 
at temperatures of 800 °C and below was pointed out 
in comparison with the annealing at 900 °C and higher 
temperatures. Really, at low annealing temperatures an 
extended "tail" is observed in the low-concentration re- 



gion of the impurity profile and the shape of this "tail" 
is a straight line if the axis of concentration is logarith- 
mic. This feature of the boron distribution is observed 
even for "tail" extension compared to the characteristic 
size of the implanted region or smaller than it. At the 
same time, after annealing at a temperature of 900 °C or 
higher the shape of the boron profile for concentrations 
below approximately 10* ^m~^ becomes convex upwards, 
i.e., similar to the Gaussian distribution. This behavior 
of the profile shape gives clear evidence of the change in 
the boron diffusion mechanism. 

It is worth noting that the "tails" that represent a 
straight line, if the axis of concentration is logarith- 
mic, are often observed directly after ion implantation of 
boron, phosphorus, gallium, and other impurity at room 
temperature (see, for example,— 1^ in the case of implan- 
tation of boron ions, ^^1^^ for phosphorus implantation, 
•^^1^^ in the case of implantation of gallium ions). The 
ion implantation in all the cases investigated was carried 
out in the direction deflecting from the crystal axis. Ac- 
cording tc4£, to completely remove the phenomenon of 
channeling and to eliminate the "tails" related to the 
scattering of ions into channels, the ion implantation 
should be carried out in amorphous silicon. Neverthe- 
less, later experiments show that "tails" are observed for 
boron implantation in the layers preamorphized by ger- 



manium ion; 
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For the low-temperature treatments of ion-implanted 
layers either an increase and broadening of existing 
"tails" in the bulk of a semiconductor occur or the for- 
mation of new "tails" if they are not observed after im- 
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plantation. In the region of low impurity concentra- 
tion, a "tail" represents as before a straight line for low- 
temperature annealing (i.e. for a small thermal budget). 
Really, an increase in the "tail" extension occurs during 
subsequent thermal treatments of preamorphized silicon 
layers that were implanted with boron ions^i^li^i^. In 
the investigations carried out b yi5i22,28,33 dearly identi- 
fied "tails" after boron implantation were not observed. 
However, such "tails" were formed in the course of the 
subsequent annealing. As follows from the experimental 
data of"*^, a "tail" is also formed in the case of thermal 
treatment at 900 °C of silicon layers implanted with in- 
dium. This "tail" represents a straight line if the concen- 
tration axis is logarithmic. The experimental data also 
show that "tails" characterized by a straight line are of- 
ten observed for the "hot" ion implantation of indium^, 
gallium^, antimony^^, and other impurities. 

It was assumed originally that the formation of "tails" 
in the low-concentration region of ion-implanted impu- 
rity profiles, especially in the case of "hot" ion implanta- 
tion, results from the fast diffusion of implanted impurity 
atoma^i^"— . It was supposed that impurity intersti- 
tials are this fast diffusing species^i^i^. For example, 
it was shown experimentally in^^ that during annealing of 
ion-implanted layers a significant fraction of indium and 
tellurium atoms leave their substitutional positions and 
become interstitials. The interstitial position is also a 
characteristic feature for atoms of gallium that, as well as 
boron, indium, and tellurium is the element of III groups. 
However, in the latest papers^i^ there are very serious 
arguments that the formation of "tails" in ion-implanted 
layers is related to the scattering of ions that reserved a 
part of the kinetic energy into channels. For example, 
in the paper— phosphorus ions were implanted into thin 
silicon layers of different thicknesses. These layers were 
located on a substrate, which collected ions channeling 
through the layer. It was assumed that diffusing atoms 
do not have sufficient energy to leave the silicon crystal, 
whereas the channeling ions have. Experiments showed 
that the substrate really collects ions passed through the 
silicon layer. The doses of the passed ions were obtained 
as a function of the layer thicknesses. These doses cor- 
respond to the doses of phosphorus atoms in the region 
of the remainder of the "tail" for the investigated depth 
if ions were implanted in the continuous silicon. The 
results obtained were generalized in^ for the cases of 
indium and gallium implantation. According to^^ the 
radiation-enhanced diffusion of gallium is impossible at 
room temperature. On the other hand, a characteristic 
"tail" is observed experimentally due to the scattering of 
ions into channels. Taking into account the possible an- 
nealing of damages and recovery of the crystal structure 
in the region of the end of the ion range, one can explain 
the results for "hot" high fluence implantation of in- 
dium and gallium ions without attracting the concept of 
an anomalous diffusion. Besides, at low temperatures a 
"tail" can be formed as a result of the channeling of a 
part of impurity atoms at the initial stage of implanta- 



tion, when the amorphous phase was not formed as yet. 
Nevertheless, the mechanism of the "tail" formation is 
not clear until now, especially taking into account the 
latest experiments22i^ related to boron implantation in 
the layer preamorphized by implantation of germanium 
ions. This allows us to formulate the following purpose 
of the research. 



II. MAIN GOAL OF THE RESEARCH 

Not rejecting the possibility for a part of the 
ions of scattering into channels, we are to show 
that the long-range migration of nonequilibrium 
impurity interstitials is the main factor in the for- 
mation of "tails" in the region of low impurity 
concentration for random ion implantation into 
silicon crystals and implantation into preamor- 
phized silicon layers. 



III. ANALYSIS OF THE MECHANISMS OF 
THE "TAIL" FORMATION DURING ION 
IMPLANTATION 

Let us consider four characteristic cases of the forma- 
tion of "tails" in the region of impurity concentration 
decreasing in the bulk of the semiconductor: 

(i) "Tail" formation during the subsequent annealing 
of ion-implanted layeva^^^^^^^^^^; 

(ii) Formation of "tails" during ion implantation into 
amorphous silicon-^'^''^l:; 

(iii) The phenomenon of "tail" formation during "hot" 
ion implantatio n^^ 

(iv) Formation of "tails" during ion implantation in 
the direction deflecting from the axis of the crystal at 
room temperature of the substrate^ '^^i^^" — . 

It is evident that in the first two cases, there is no 
phenomenon of channeling and it is possible to explain 
the formation of "tails" only by the anomalous impurity 
diffusion. It is worth noting that impurity diffusion can 
occur in an amorphous phase too, especially, if the high- 
current implanter is used for ion implantation and there 
is a possibility of substrate heating, as it was the case 
in the experiments of Cristiano et al;^ As an example. 
Fig. [T] presents the calculation of ion-implanted indium 
redistribution. The indium concentration profile was cal- 
culated within the framework of the model for diffusion 
of impurity interstitials described below. For compari- 
son, the experimental data of*^ are used. In the work of 
Gamo et alM- the distributions of impurity atoms were 
obtained by measuring the 7-ray intensities in a com- 
bination with the layer removal technique. Indium was 
implanted with an energy of 45 keV to a dose of ~1 x 10^^ 
ion/cm^ and 8° off the < III > axis in order to reduce 
channeling effects. The temperature of annealing was 900 
°C, and the thermal treatment duration was 20 minutes. 
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The following values of parameters for the model of in- 
terstitial diffusion of ion-implanted impurity were used. 
The parameters of the as-implanted indium distribution 
are: Rp = 0.03 fim (30 nm); ARp = 0.0084 ^m (8.4 nm); 
the parameters of the indium interstitial diffusion are: 
the average migration length of indium interstitials Iai 
— 0.042 /im (42 nm); the time-average value of the gener- 
ation rate of nonequilibrium impurity interstitials in the 
maximum of distribution = 7.9x10^ fim~^s~^. Here 
Rp and ARp are the average projective range of impurity 
ions and straggling of the projective range, respectively. 
These parameters provide the best fitting of the calcu- 
lated indium concentration profile to the experimental 
one. 

As can be seen from Fig. [1] the results of calculation 
agree well with the experimental data. It is a very impor- 
tant argument in favour of the fast diffusion of nonequi- 
librium indium interstitials. 
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FIG. 1: Calculated indium concentration profile (solid line) 
after thermal treatment of implanted silicon substrate at 900 
°C for 20 minutes. The experimental data (x - impurity 
distribution after implantation; • - after annealing) are taken 
frorr4i and the indium solubility limit is taken from — 



Taking into account the good agreement with the ex- 
perimental data, we use the model of the interstitial dif- 
fusion of ion-implanted impurity for simulating indium 
redistribution during "hot" ion implantation that was in- 
vestigated in"*^ . The results of simulation are presented 
in Fig. [2j The energy and the dose of indium ions are 
the same as in the previous experiment, namely, 45 keV 
and ^1x10^^ ions/cm'^. The following value of the aver- 
age migration length of indium interstitials was used for 
the best fit to the experimental data: Iai — 0.032 /im 
(32 nm). This value is greater than the average migra- 
tion length Iai — 0.027 fj,m used for similar calculations 
in^. Perhaps, the increase of Iai is due to the effect 
of evaporation of indium interstitials from the surface of 
the semiconductor^^ that we took into consideration in 



contrast to^^. Taking into account the evaporation of 
indium interstitials made it possible to achieve the best 
agreement with the experimental profile in the vicinity 
of the surface. 
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FIG. 2: Calculated indium concentration profile (solid line) 
after "hot" ion implantation. The experimental data are 
taken from^ 



As can be seen from Fig. [2J the results of calculation 
agree well with experimental data. The value of the av- 
erage migration length for "hot" implantation of indium 
ions Iai = 0.032 fim is close to the migration length Iai 
= 0.042 /im for annealing of indium implanted layer at a 
temperature of 900 °C. It means that in the experiment 
with "hot" ion implantation a "tail" is also formed due 
to the fast migration of impurity interstitials. 

The results obtained suggested to revise the inferences 
drawn in^i^ about "tail" formation due to the scattering 
of ions into channels. Really, these inferences were based 
on the experiments22 with phosphorus ion implantation 
in thin silicon layers of different thicknesses. The main 
assumption of the above-mentioned paper is that diffus- 
ing phosphorus atom cannot pass through the boundary 
of the silicon crystal, while the ions which scatter into 
channels can cross the boundary due to the rest of the 
kinetic energy. However, within the framework of the in- 
terstitial diffusion model the "tail" formation is related 
to the diffusion of impurity interstitials. Therefore, there 
are no restrictions for impurity interstitials to evaporate 
from the surface of a crystal and then to collect on the 
substrate, on which the implanted layer is placed. This is 
confirmed by the experiments concerning hydrogen diffu- 
sion through thin silicon layers^i, because the extended 
"tails" on the hydrogen concentration profiles are also 
formed by interstitial diffusion^. 

It is clear from the analysis carried out above that 
the vast majority of phenomena of "tail" formation in 
the low-concentration region of ion-implanted impurity 
profile is due to the fast impurity diffusion. The scat- 
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tering of ions into channels can play only an auxiliary 
role in the experiments under investigation. It is usually 
supposed that impurity interstitials are the fast diffusing 
species^i^i^'^. Let us note that the mechanism of the 
long-range migration of nonequilibrium impurity intersti- 
tials was already used for explanation of gold diffusion in 
ion-implanted layersfiS. This mechanism was also used for 
describing "tail" formation during thermal treatments of 
the layers implanted with boro n^^'^^ and for explanation 
of boron diffusion in "buried" layers under the conditions 
of generation of nonequilibrium silicon interstitials^^iS. 
Really, in^ an analytical solution describing the redis- 
tribution of impurity atoms due to diffusion of boron in- 
terstitials during thermal treatment was obtained. As 
the initial condition for the distribution of boron atoms 
the (5-function was chosen. It was shown that for single 
participation of a boron atom in the act of interstitial 
migration a "tail" represents a straight line if the axis 
of concentration is logarithmic. It is necessary to note 
that for the initial boron distribution in the form of S- 
function, the requirement of the occurrence of one event 
of boron interstitial migration is equivalent to the con- 
dition of the long-range migration of boron interstitials. 
On the other hand, under the condition of multiple mi- 
grations of the same impurity atom the boron profile has 
the form of Gaussian distribution^^,. Therefore, it was 
proposed in'^'^ to determine the microscopic mechanism 
of impurity transport by means of a change in the form 
of the impurity profile with increase in annealing dura- 
tion. It is worth noting that the calculations which were 
carried out 35^48,49^54 gj^Q^ ^]^g^^ ^j^g "tail" represent- 
ing a straight line is formed by the long-range migration 
of impurity interstitials not only in the case of initial dis- 
tribution in the form of (5-function, but also for a genera- 
tion rate of nonequilibrium interstitials described by the 
Gaussian distribution or any other distribution. Besides, 
the average migration length of impurity interstitials can 
be comparable with the characteristic size of the initial 
doped region, or even less than this size. 



us consider for definiteness that boron transport during 
low-temperature annealing is carried out by migration of 
nonequilibrium impurity interstitials. We shall neglect 
the diffusion of the "impurity atom - intrinsic point de- 
fect" pairs in view of the small thermal budget. This 
assumption is based on the well-known opinion that im- 
purity interstitials are the most mobile species. Taking 
into account that the extended "tail" is formed in the 
low-concentration region, where C < , we can suppose 
that the electric field does not infiuence the migration of 
impurity interstitials. Here C is the concentration of the 
substitutionally dissolved boron atoms; Ui is the intrin- 
sic carrier concentration for the treatment temperature. 
It is clear that there is also no infiuence of the electric 
field in the case of migration of neutral interstitials even 
if C > Tii. We shall also neglect the infiuence of elas- 
tic stresses and suppose that substitutionally dissolved 
boron atoms and boron atoms incorporated into clusters 
and extended defects are immobile. Then, the system 
of equations describing redistribution of ion-implanted 
boron during low-temperature annealing has the follow- 
ing form: 

(i). The conservation law for immobile impurity 
atoms: 



dt 



rAI 



(1) 



(ii). The equation of diffusion for nonequilib- 
rium impurity interstitials: 



q2 /^AI r<AI 



or 



92 c 



AI 



c 



AI 



'A/ 



~g^'{x,t) 



'A/ 



(2) 



(3) 



where 



IV. MODEL OF INTERSTITIAL DIFFUSION 

In the paper^^, an equation for migration of nonequi- 
librium impurity interstitials was obtained which takes 
into account the different charge states of migrating 
species and also the drift of impurity interstitials in the 
built-in electric field and in the field of the elastic stresses. 
This equation has a form similar to the form of the equa- 
tion describing the diffusion of nonequilibrium "impurity 
atom - intrinsic point defect" pairs when the influence 
of the built-in electric field and the field of the elastic 
stresses are also taken into account. As the mathematical 
formulations for the migration of nonequilibrium impu- 
rity interstitials and for the diffusion due to the nonequi- 
librium pairs are identical, it is impossible to make an ex- 
act conclusion about the kind of the mobile species from 
the form of impurity profile after thermal treatment. Let 



Iai = VdAh^, 



~g^\x,t)^G^\x,t)T 



AI 



(4) 



Here is the total concentration of substitutionally 
dissolved impurity atoms and impurity atoms incorpo- 
rated into clusters or trapped by extended defects (im- 
mobile impurity atoms) : C^^ is the total concentration of 
nonequilibrium impurity interstitials in different charge 
states; d-^^ and t^^ are the diffusivity and the average 
lifetime of nonequilibrium impurity interstitials, respec- 
tively; G^^ is the generation rate of impurity intersti- 
tials. We use the stationary diffusion equation for im- 
purity interstitials in view of their large migration length 
Iai (Iai ^ I fail , where I faii is the characteristic length of 
the decrease in the impurity concentration in the high- 
concentration region of impurity profile), and due to a 
small lifetime of these nonequilibrium interstitials r'^^ 
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(r <^ Tp, where Tp is the duration of anneahng). Let 
us note that the system of equations ([T]) , ([2]) is similar to 
the systems of equations used in^"— for the description 
of interstitial diffusion. 

To describe the spatial distribution of impurity atoms 
after implantation and spatial distribution of the genera- 
tion rate of boron intcrstitials, the Gaussian distribution 
is chosen: 



Co (a;) = C{x, 0) = C,„ exp 



{x - Rp)^ 
2AR^ 



G^^(a:,0 =5™ exp 



(x - Rp) 



21 



where 



Q 



X 10" 



2Ai?2 



[/im 



(5) 



(6) 



(7) 



Here Cm is the maximal concentration of impurity 
atoms after implantation; is the maximal value of 
generation rate of impurity intcrstitials per unit volume; 
Q is the dose of ion implantation [ion/cm^]; Rp and Ai?p 
are the average projective range of impurity ions and 
straggling of the projective range, respectively. 

It is clear from expressions ([5]) and ^ that the model 
presupposes that the generation rate of boron intcrsti- 
tials is proportional to the total concentration of impurity 
atoms. Really, the analysis carried out ini^ shows that 
the boron intcrstitials can be generated during thermal 
treatment not only as a result of annealing of radiation 
defects, but also due to dissolution or rearrangement of 
the clusters that incorporated impurity atoms and also as 
a result of elastic stresses arising due to the small atomic 
radius of boron. This assumption is confirmed by the ex- 
perimental data of^, where the long-range migration of 
boron interstitials with the formation of the characteris- 
tic "tail" was observed during doping of the amorphous 
silicon layer deposited with the low thermal budget on a 
silicon substrate and then subjected to solid phase epi- 
taxy. Thus, it follows front^ that the long-range migra- 
tion of boron interstitials is observed in the case of no 
radiation defects. As both the cluster concentration and 
intensity of elastic stresses are proportional to the total 
concentration of impurity atoms, the use of expression (6) 
for describing the generation rate distribution of boron 
interstitials is quite reasonable. 



V. RESULTS OF SIMULATION OF 
ION-IMPLANTED BORON REDISTRIBUTION 

In the case of the constant coefhcients and the system 
of equations ([l} and Q allows obtaining an analytical so- 
lution. Such solutions for different boundary conditions 



on a finite interval [0,Xi?] within the assumption of con- 
tinuous generation of impurity interstitials in the doped 
layer were obtained in^i^. We used these solutions for 
modeling the redistribution of ion-implanted boron. Let 
us note that the diffusivity d^^ has a constant value, if 
the migration of boron interstitials is not influenced by 
the built-in electric field. The average lifetime of these 
interstitials has a constant value in the case of absorption 
by the unsaturated traps distributed homogeneously in 
the bulk of a semiconductor. For example, such case of 
absorption is realized if replacement of the silicon atom 
by the boron interstitial from the lattice site to the in- 
terstitial one occurs. 

Let us consider in the beginning the results of simula- 
tion of the experimental data o^, because in this inves- 
tigation various methods of the transient enhanced dif- 
fusion suppression were used. So, in^ for suppression of 
the transient enhanced diffusion Czochralski grown (100) 
n-type silicon wafers were subjected to preamorphization 
by performing a germanium (Ge) ion implantation at 15 
keV (Ge-PAI) and subsequently with 1 keV boron ions to 
a dose of 1.5x10^^ cm"^. Nitrogen (N) co-implantation 
to the same dose of 1.5x10^^ cm^'^ was performed on 
some wafers. The thermal annealing was carried out in 
a rapid thermal processing system under N2 ambient du- 
ration at a temperature of 800 °C for 60 seconds. The 
dopant profiles were analyzed ex-situ by secondary ion 
mass spectrometry (SIMS). A primary Cross-sectional 
TEM (XTEM) was also performed to analyze the extent 
of amorphization and the Ge-PAI induced EOR defects. 
The simulation of boron redistribution carried out for 
these conditions in'^^, has shown good agreement with 
the experimental data for the average migration length 
of boron interstitials Iai = 11 nm. Also it was sup- 
posed that approximately 8.6% of the implanted boron 
atoms occupied interstitial positions. Migration of these 
nonequilibrium interstitial atoms results in the formation 
of an extended "tail" on the boron concentration profile. 
This "tail" is located in the interval from 0.02 /im up 
to approximately 0.1 /im from the surface of a semicon- 
ductor. For comparison Fig. [3] presents the results of 
modeling for the same process of ion-implanted boron 
redistribution in the case of additional nitrogen implan- 
tation at 6 keV after preamorphization due to germanium 
ions'^'^. The energy of nitrogen implantation equal to 6 
keV was chosen so that the maximum of N concentration 
was located between the boron doped region and a/c in- 
terface. 

As can be seen from Fig. [3l there is good agreement 
between the calculated profile and the experimental data. 
The following values of the model parameters were used 
to provide the best fit of the calculated boron concen- 
tration profile to the experimental one: Parameters 
prescribing the initial distribution of implanted 
boron: Rp = 0.0044 //m (4.4 nm); Ai?p = 0.0036 /im 
(3.6 nm). Parameters specifying the process of in- 
terstitial diffusion: the duration of annealing Tp — 60 s; 
the annealing temperature T = 800 °C; the average life- 
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FIG. 3: Calculated boron concentration profile in silicon 
preamorphized by Ge ions after annealing for 60 s at a temper- 
ature of 800 °C. The dotted curve represents calculated boron 
distribution after ion implantation at an energy of 1 keV. An 
additional implantation of nitrogen ions at an energy of 6 
keV for stronger suppression of the transient enhanced diffu- 
sion was carried out. The experimental data (filled circles) 
are taken from the papep^ 



time of nonequilibrium boron interstitials t^^ = 0.011 s; 
the maximum value of the generation rate of nonequilib- 
rium impurity interstitials = 6.7 x 10^ fim.~^s~^; 
the average migration length of boron interstitials Iai 
— 5.9 nm; the concentration of boron interstitials on the 
right boundary Cp^ — 0; the position of the right bound- 
ary xp = 0.5 fim. It was supposed that approximately 21 
% of the implanted boron atoms occupied the interstitial 
positions temporally. 

As can be seen from the calculations performed, the 
additional implantation of nitrogen ions results in a sig- 
nificant (almost 2 times) reducing of the average migra- 
tion length of boron interstitials. Therefore, the shrink- 
age of the "tail" approximately by 0.03 fim is observed 
that is very attractive from the technological point of 
view. It is worth noting that there is a significant (2.8 
times) increase in the generation rate of boron intersti- 
tials. However, due to the reduced migration length of 
these interstitials, their main fraction is trapped within 
the implanted layer that does not result in an increase of 
the depth of p — ti junction. We assume that this trap- 
ping is due to the interaction with nitrogen atoms or with 
complexes of nitrogen and germanium atoms. 

A similar modeling for the case of transient enhanced 
diffusion suppression by using the lower annealing tem- 
perature, namely 750 °C, is presented in Fig. 01 The 
following values of the parameters specifying the 
process of interstitial diffusion were used: the max- 
imum value of the generation rate of nonequilibrium im- 
purity interstitials g:^^ = 7.4 x 10^ ^m^^s^^; the aver- 
age migration length of boron interstitials Iai — 4.7 nm. 



It follows from the values obtained that the decrease of 
the annealing temperature results in the reduction of the 
average migration length of boron interstitials and in the 
highly abrupt impurity profile. 



c 

(U 

o 
c 
o 
o 

c 

e 

o 

"to 
o 



10 

10' 

10^ 

10' 

10' 

10^ 




— calculated 

• experiment 

• - - initial 

— solubility 



— intrinsic carrier 
concentration 



_L 



J 



0.00 0.02 0.04 0.06 0.08 0.10 
Depth (^m) 



FIG. 4: Calculated boron concentration profile in silicon 
preamorphized by Ge ions after annealing for 60 s at a temper- 
ature of 750 °C. The dotted curve represents calculated boron 
distribution after ion implantation at an energy of 1 keV. An 
additional implantation of nitrogen ions at an energy of 6 
keV for stronger suppression of the transient enhanced diffu- 
sion was carried out. The experimental data (filled circles) 
are taken from the papep^ 



Now, let us present the results of modeling the in- 
terstitial diffusion during boron implantation using the 
data obtained in^^. In the paper~, n-type 10-20 flcm, 
prime (100) Si wafers were used. The wafers were first 
preamorphized with Ge ions (30 keV, 1 x 10^^ cm~^) 
to a depth of about 50 nm and then implanted with 0.5 
keV boron ions to a dose of 1 x 10^^ cm~^. The im- 
plants were performed on a high-current implanter at tilt 
and twist angles of 0° with electrostatic deceleration in 
front of the target. It is worth noting that due to the 
high-current implanter the wafer can be heated during 
ion implantation. The dopant atom distributions of all 
samples were analyzed using a high-resolution secondary 
ion mass spectrometry (SIMS). The boron concentration 
profile measured after implantation in preamorphized sil- 
icon is shown in Fig. [5l 

As can be seen from Fig. [51 the calculated curve agrees 
well with the measured boron concentration profile after 
implantation characterized by the extended "tail" . Thus, 
the experimental data^'^ can be explained on the basis of 
the long-range migration of nonequilibrium boron atoms. 
Really, there is no channeling of ions due to implantation 
in the amorphous phase, whereas fast anomalous diffu- 
sion is quite possible in amorphous silicon. The following 
values of the model parameters were used to provide the 
best fit of the calculated boron concentration profile to 
the experimental one: Parameters prescribing the 
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FIG. 5: Calculated concentration profile of boron implanted 
into silicon preamorphized by Ge ions. The long-range mi- 
gration of boron atoms during implantation is taken into ac- 
count. The dotted curve represents boron distribution after 
ion implantation at an energy of 0.5 keV calculated with no 
interstitial diffusion. The experimental data (filled circles) are 
taken from^S 



0.0014 fim (1.4 nm). Parameters specifying the pro- 
cess of interstitial diffusion: the maximum value of 
the time-average generation rate of nonequilibrium impu- 
rity interstitials g^/ — 2.0x10^ fim^^s^^; the average 
migration length of boron interstitials Iai = 0.0072 /im 
(7.2 nm). It was supposed that approximately 4.1 % of 
the implanted boron atoms occupied the interstitial po- 
sitions temporally. 



10' 



c 
o 



10^ 



r 10" r- 



C 
(U 

o 
c 
o 
o 

c 
o 

o 
"to 

-I— » 

o 



— calculated 

• experiment 
■ ■ - initial 

• — solubility 



10' 



10° 



10 




— intrinsic carrier 
concentration 



0.00 



0.02 0.04 
Deptli (nm) 



0.06 



Gaussian distribution of implanted boron: Rp = 

0.0006 fim (0.6 nm); Ai?p = 0.0027 (2.7 nm). Pa- 
rameters specifying the process of fast diffusion: 

the maximum value of the generation rate of nonequi- 
librium impurity interstitials g:^^ — 3.2x10^ ^m~'^s~^; 
the average migration length of boron interstitials Iai = 
0.01 iim (10 nm). For the best fit it was supposed that 
approximately 0.63 % of the implanted boron atoms par- 
ticipated in the long-range migration. However, there is 
also another way to explain the "tail" formation'^''. Re- 
ally, the boron ions were implanted with electrostatic de- 
celeration in front of the target. During the deceleration 
stage, some ions can become neutral due to reciprocal 
collisions. In such a case they will not be decelerated 
and keep their original energy. Due to the higher energy, 
these ions penetrate deeper in the bulk of the semicon- 
ductor. Therefore, this problem requires a further inves- 
tigation. 

Finally, in Fig. |6] the results of modeling of ion- 
implanted boron redistribution are shown, when the im- 
plantation of BF2 ions was used for amorphization of the 
surface layer of silicon crystal and transient enhanced 
diffusion suppression^. In^, (100) oriented, 10 ilcm, n- 
type silicon wafers were subjected to BF2 implantation 
at 2.2 keV to a dose of l.Ox 10^^ cm^^. The furnace an- 
nealing was carried out at a temperature of 800 °C for 30 
minutes. The resulting boron concentration profiles were 
measured by SIMS. The following values of the model 
parameters were used to provide the best fit of the calcu- 
lated boron concentration profile to the experimental one: 
Parameters prescribing the initial distribution of 
implanted boron: Rp — 0.0002 fim (0.2 nm); ARp = 



FIG. 6: Calculated boron concentration profile after anneal- 
ing of a silicon layer implanted with BF2 ions for TED sup- 
pression. The dotted curve represents the calculated boron 
distribution after implantation of BF2 ions. The annealing 
was carried out for 30 minutes at a temperature of 800 °C. The 
experimental data (open circles after implantation and filled 
circles after annealing) are taken from the paper by Suzuki^ 

As is seen from Fig. [HI the calculated boron profile 
after furnace annealing agrees well with the experimen- 
tal one in the "tail" region of impurity distribution that 
look like a straight line. It means that the "tail" region 
is formed due to the long-range migration of nonequilib- 
rium boron interstitials. It is necessary to note that mod- 
eling of the experimental data obtained in^ was carried 
out earlier in the paper- ^ . The model of Uematsi*^ was 
used for simulation of boron diffusion ir*^. It is supposed 
in this model that there is a local thermodynamic equi- 
librium between substitutionally dissolved boron atoms, 
silicon self-interstitials, and boron interstitials. It means 
from the microscopic point of view that the same boron 
atom becomes frequently interstitial, makes a number of 
chaotic jumps, and then becomes substitutionally dis- 
solved again due to the "kick-out" of silicon atom from 
the lattice site. According tc^^, the impurity profile for 
this diffusion mechanism has a form of the Gaussian dis- 
tribution. Really, the boron concentration profile calcu- 
lated irt^ represents the Gaussian distribution and dis- 
agrees with the experimental data. On the other hand, 
as follows from the experimental data^*, the boron con- 
centration profile after annealing for 1.5 h begins to get 
a convex form similar to the Gaussian distribution. Re- 
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maining within the framework of the model of the long- 
range migration, such a change in the form of impu- 
rity distribution can be explained by taking into account 
the significant increase in the boron concentration in the 
"tail" region. Really, if the impurity concentration is 
greater than the intrinsic carrier concentration n^, the mi- 
gration of charged boron interstitials is influenced by the 
built-in electric field. On the other hand, the increase of 
annealing duration can create appropriate conditions for 
intensive dissolution of the clusters and extended defects 
incorporating silicon self-interstitials. In turn the gener- 
ation of nonequilibrium self-interstitials results in boron 
diffusion due to the generation, migration, and dissolu- 
tion of the "boron atom - silicon self-interstitial" pairs. 
It is usually supposed (see, for example,—"—) that the lo- 
cal thermodynamic equilibrium prevails between substi- 
tutionally dissolved boron atoms, silicon self-interstitials, 
and the pairs. Therefore, the boron concentration pro- 
file formed due to the pair diffusion should look similar 
to the Gaussian distribution that agrees with the experi- 
mental data for large thermal budget. However, the final 
solution of this problem requires further careful investi- 
gations. 

VI. CONCLUSIONS 

We have shown that the vast majority of the phenom- 
ena of "tail" formation in the low-concentration region 
of ion-implanted impurity profile, including ion implan- 



tation in the preamorphized silicon, "hot" ion implanta- 
tion, impurity redistribution during the subsequent low- 
temperature annealing of ion-implanted layers, is related 
to the fast impurity diffusion. The channeling of ions 
scattered into channels is negligible for all the cases under 
consideration. With this purpose, simulation of indium 
redistribution during the subsequent annealing of ion- 
implanted layers and during "hot" indium implantation 
has been carried out. Many cases of the "tail" formation 
under different conditions of transient enhanced diffusion 
suppression during low-temperature annealing of boron 
implanted layers have been also investigated, including 
boron implantation into silicon preamorphized by ger- 
manium ions. The impurity concentration profiles calcu- 
lated within the framework of the model of the long-range 
migration of nonequilibrium impurity interstitials agree 
well with experimental data. On this basis the values of 
the average migration length of nonequilibrium impurity 
interstitials have been obtained. So, if for boron implan- 
tation in silicon crystals preamorphized by germanium 
ions, the average migration length of boron interstitials 
is equal to 11 nm at an annealing temperature of 800 °C, 
the additional nitrogen implantation allows one to reduce 
this value approximately to 6 nm. A further reduction of 
the average migration length of boron interstitials can be 
achieved by decreasing the annealing temperature to 750 
°C. Finally, for BF2 implantation in silicon crystals the 
average migration length of boron interstitials is equal to 
7.2 nm at an annealing temperature of 800 °C. 
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